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We analyze the perturbative and parametric stability of the QCD predictions for the Callan-Gross 
ratio R{x,Q^) = Fl/Fj and azimuthal cos,{2(p) asymmetry in heavy-quark leptoproduction. Our 
analysis shows that large radiative corrections to the structure functions cancel each other in their 
ratio R[x^ Q^) and azimuthal asymmetry with good accuracy. As a result, the NLO contributions to 
the Callan-Gross ratio and cos(2(p) asymmetry are less than 10% in a wide region of the variables 
X and 2^. We provide compact analytic predictions for R{x,Q^) and asymmetry in the case of 
low X <C 1 . Simple formulae connecting the high-energy behavior of the Callan-Gross ratio and 
azimuthal asymmetry with the low-x asymptotics of the gluon density in the target are derived. 
It is shown that the obtained hadron-level predictions for R{x, Q^) and azimuthal asymmetry are 
stable at X ^ 1 under the DGLAP evolution of the gluon distribution function. 
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1. Introduction 

In the framework of perturbative quantum chromodynamics (QCD), the basic spin-averaged 
characteristics of heavy-flavor photo- [|l], electro- [||], and hadro-production [^-^ are known 
exactly up to the next-to-leading order (NLO).^ Although these explicit results are widely used at 
present for a phenomenological description of available data (for a review, see Ref. [||]), the key 
question remains open: How to test the applicability of QCD at fixed order to heavy-quark pro- 
duction? The basic theoretical problem is that the NLO corrections are sizeable; they increase the 
leading-order (LO) predictions for both charm and bottom production cross sections by approxi- 
mately a factor of two. Moreover, soft-gluon resummation of the threshold Sudakov logarithms 
indicates that higher-order contributions can also be substantial. (For reviews, see Refs. [ |l0| , [TTI].) 
On the other hand, perturbative instability leads to a high sensitivity of the theoretical calculations 
to standard uncertainties in the input QCD parameters. The total uncertainties associated with the 
unknown values of the heavy-quark mass, m, the factorization and renormalization scales, pLp and 
}Xr, the asymptotic scale parameter Aqcd and the parton distribution functions (PDFs) ai^e so large 
that one can only estimate the order of magnitude of the pQCD predictions for cham production 



cross sections in the entire energy range from the fixed-target experiments [|TJ, 13 1 to the RHIC 
colhder [|]. 

Since these production cross sections are not perturbative ly stable, it is of special interest to 
study those observables that are well-defined in pQCD. Nontrivial examples of such observables 



were proposed in Refs. [14-22], where the azimuthal cos (29) asymmetry and Callan-Gross ratio 
R{X)Q^) = Fl/Pt in heavy-quark leptoproduction were analyzed.^'^ In particular, the NLO soft- 
gluon corrections to the basic mechanism, photon-gluon fusion (GF), were calculated. It was shown 
that, contrary to the production cross sections, the azimuthal cos(2<p) asymmetry in heavy-flavor 
photo- and leptoproduction is quantitatively well defined in pQCD: the contribution of the domi- 
nant GF mechanism to the asymmetry is stable, both parametrically and perturbatively. Therefore, 
measurements of this asymmetry should provide a clean test of pQCD. 

The perturbative and parametric stability of the GF predictions for the Callan-Gross ratio 
R{xiQ^) = Fl/Fj in heavy-quark leptoproduction was considered in Refs. [|o|-^. It was shown 
that large radiative corrections to the structure functions Ft{x,Q^) and Fl{x,Q^) cancel each other 
in their ratio R{x,Q^) with good accuracy. As a result, the next-to-leading order (NLO) contribu- 
tions of the dominant GF mechanism to the Callan-Gross ratio are less than 10% in a wide region 
of the variables x and Q^. 

In the present paper, we continue the studies of perturbatively stable observables in heavy- 
quark leptoproduction, 

+N{p) ^ l{l -q) + Q{pq) +X[Q]{px). (1.1) 



' Some recent results concerning the ongoing computations of tire next-to-next-to-leading order (NNLO) corrections 
to the heavy-flavor hadroproduction are presented in Refs. ^ 

^Well-known examples include the shapes of differential cross sections of heavy flavor production, which are suffi- 
ciently stable under radiative corrections. 

■'Note also the paper |^3|], where the perturbative stability of the QCD predictions for the charge asymmetry in 
top-quark hadroproduction has been observed. 
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Figure 1: Definition of the azimuthal angle (p in the nucleon rest frame. 



In the case of unpolaiized initial states and neglecting the contribution of Z-boson exchange, the 



azimuth-dependent cross section of the reaction ( |1.1| ) can be written as 



2< 



Ft {x, ) + eFiXx, Q^) + eFa 2^ ) cos 2<p + 2 ^/e{l + e)Fi{x, ) cos <p 

(1.2) 

where aem is Sommerfeld's fine-structure constant, F2{x,Q^) = 2x{Ft +Fl), the quantity e mea- 
sures the degree of the longitudinal polarization of the virtual photon in the Breit frame [El], 



2(1 -v) 
l+{l->t 



, and the kinematic variables are defined by 



p-q 



2p-q' 



xyS, 



m 



2 • 



(1.3) 



In Eq. (1.2), Ft{Fi) is the usual "fN structure function describing heavy-quark production by 
a transverse (longitudinal) virtual photon. The third structure function. Fa, comes about from 
interference between transverse states and is responsible for the cos2<p asymmetry which occurs in 
real photoproduction using linearly polarized photons. The fourth structure function, Fj, originates 



from interference between longitudinal and transverse components []24[]. In the nucleon rest frame, 
the azimuth (p is the angle between the lepton scattering plane and the heavy quark production 
plane, defined by the exchanged photon and the detected quark Q (see Fig. |I|). The covariant 
definition of <p is 



cos 9 



r-n 



Pv^a^Pj 



Q^Jl/x^+4ml/Q^ 

sirup = — — n-i 

2v — r^V — « 

n^=e^'''Pq,PaPQp. 



(1.4) 
(1.5) 



In Eqs. ( |1.3[ ) and ( |1.5[ ), m and are the masses of the heavy quark and the target, respectively. 

In this talk, we review the perturbative and parametric stability of the Callan-Gross ratio, 
R{x,Q^), and azimuthal cos(2(p) asymmetry, A{x,Q^), defined as 



R{x,Q^) = ^{x,Q^) 
rj 



A{x,Q' 



2x^ix,Q^) 



(1.6) 
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First, we consider radiative corrections to the quantity R{x,Q^) using the explicit NLO results 
presented in 25|. Our calculations show that complete ^{ocj) corrections to R{x, Q^) (including 
both the photon-gluon, y*g — )• QQ.{g), and photon-(anti)quark, Y'q — )• QQq, fusion components) do 
not exceed 10% in the energy range x > 10^^. 

Then, we analyze the perturbative stability of the azimuthal cos(2<p) asymmetry, A{x,Q^). 
Presently, the exact NLO predictions for the azimuth dependent structure function Fa{x, Q^) are not 
available. For this reason, we use the so-called soft-gluon approximation to estimate the radiative 
corrections to Fa{x, Q^). Our analysis shows that the NLO soft-gluon predictions for A(;c, Q^) affect 
the LO results by less than a few percent at < and x > 10^^. 

In both cases, perturbative stability is mainly due to the cancellation of large radiative correc- 
tions to the structure functions Fi, Ft, Fa and F2 in their ratios, R{x, Q^) and A{x, Q^), con^espond- 
ingly. Note also that both the LO and NLO predictions for the Callan-Gross ratio and azimuthal 
asymmetry are sufficiently insensitive, to within ten percent, to standard uncertainties in the QCD 
input parameters jXp, jiR, Aqcd and PDFs. 

We conclude that, in contrast to the production cross sections, the ratios R{x, Q^) and A(;c, Q^) 
in heavy-quark leptoproduction are observables quantitatively well defined in pQCD. Measure- 
ments of these quantities in charm and bottom leptoproduction should provide a good test of the 
conventional parton model based on pQCD. 

Since the ratios R{x^Q^) and A{x,Q^) are perturbatively stable, it makes sense to provide 
the LO hadron-level predictions for these quantities in analytic form that may be useful in some 
applications. For this reason, we derive compact hadron-level LO predictions for the the Callan- 
Gross ratio and azimuthal asymmetry in the limit of low x ^ 0. Assuming the low-x asymptotic 
behavior of the gluon PDF to be of the type g{x, Q^) oc 1 /jc'+^, we provide analytic result for the 
ratios R{x — 0, Q^) and — 0, Q^) for arbitrary values of the parameter 5 in terms of the Gauss 
hypergeometric function.^ 

In principle, the parameter 5 is a function of and this dependence is calculated using the 
DGLAP evolution equations [^-Q. However, our analysis shows that hadron-level predictions 
for R{x ^ 0,0^) and A{x ^ 0,Q^) are practically independent of d in the entire region of for 
5 > 0.2. We see that the hadron-level predictions for R{x — )■ 0, 2^) and — > 0, Q^) are stable not 
only under the NLO corrections to the partonic cross sections, but also under the DGLAP evolution 
of the gluon PDF. 

As to the experimental applications, we show that our compact LO formulae for R{x — 0, Q^) 
conveniently reproduce the HERA results for fj^jc, 2^) and fj'(jc,2^) obtained by HI Collab- 
oration [33, 34] with the help of more cumbersome NLO estimations of Fl{x,Q^). Our analytic 
predictions will also be useful in extraction of the azimuthal asymmetries from the incoming COM- 
PASS results as well as from future data on heavy-quark leptoproduction at the proposed EIC [35] 
and LHeC [36] colliders at BNL/JLab and CERN, correspondingly. 

This paper is organized as follows. In Section ^ we analyze the exact NLO results for the 
Callan-Gross ratio. The soft-gluon contributions to A{x,Q^) are investigated in Section ^. The 
analytic LO results for the ratios R{x, Q^) and Q^) at low x are discussed in Section 0. 



^The simplest case, 5 = 0, has been studied in Ref. [Q. The choice 5=1/2 historically originates from the BFKL 
resummation of the leading powers of In ( 1 /x) [p7|-p^. 
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2. Exact NLO predictions for the Callan-Gross ratio R{x, Q^) 

At leading order, ^(aemOti)> leptoproduction of heavy flavors proceeds through the photon- 
gluon fusion (GF) mechanism, 



f{q)+g{kg)^Q{pQ)+Q{pQ). 



(2.1) 



The relevant Feynman diagrams are depicted in Fig. The corresponding y*8 cross sections, 

Kg 



af Vz, A) {k = 2,L,A,I), have the form 



4%-^) = g^«(z){[(l-z)^ + z2 + 4Az(l-3z)-8AV]ln|±| 
-[l+4z(l-z)(A-2)]^,}, 





X) = 


K ^ 




X) = 


-as(z)z|2A 




X) = 


0, 


= {2n] 


^ggaemz/G^, where 



(2.2) 



charge and a., = 05i(/x|) is the strong-coupling constant. In Eqs. (22), we use the following defini- 
tion of partonic kinematic variables: 



2q-kg 



X 



m 



AXz 



(2.3) 



The hadron-level cross sections, Ok,GF{x,Q^) (k = 2,L,A,/), coiTesponding to the GF subprocess, 
have the fomi 



x-(l+4A) 

where g{z,IJ-F) is the gluon PDF of the proton. 



Ok,GF{x,Q^) = I dzg{z,lJ-F)dk,g{x/z,X,iXF): 

Ja-(1+4A) 



(2.4) 





Figure 2: Feynman diagrams of photon-gluon fusion at LO. 
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The leptoproduction cross sections Ok{x, Q^) are related to the structure functions Q^) as 
follows: 



Fkix^Q") = ^ f a,(^,(2') {k = T,L,A,I), 

~t2 



F2{x,Q^) = t3^<^2{x,Q^). (2.5) 

471 Mem 

where a2(x, Q^) = Ot{x, Q^) + Ol{x, Q^). 

At NLO, i^(aem«5^), the contributions of both the photon-gluon, 'fg — > QQ{g), and photon- 
(anti)quark, Y'q — QQq, fusion components are usually presented in terms of the dimensionless 
coefficient functions c["''^(z,A) as 



nr ' 



(1,0)/ 1 ^ , (1,1)/ ^ X , M 

4 ^(z,A)+4^(z,A)ln-2 



+ ^(a2) , (2.6) 



where we identify ^ = jif = 

The coefficients c[Y^(z, A) and c[Y^(z,A) {k = T,L,A,I) of the -dependent logarithms can 
be evaluated explicitly using renomialization group arguments ^ |3|]. The results of direct calcula- 
tions of the coefficient functions c[^„°^(z, A) and cS/'"''(z, A) {k = T,L) are presented in Refs. [S, El]. 



k,g V<o'V '-/t,^ 

Using these NLO predictions, we analyze the dependence of the ratio R{x, Q^) = Fi^/Fj at fixed 
values of x. 

The panels {a), (b) and (c) of Fig. || show the NLO predictions for Callan-Gross ratio R{x, Q^) 
in charm leptoproduction as a function of ^ = Q^/m^ atx= 10^', 10^^ and 10^^, correspondingly. 
In our calculations, we use the CTEQ6M parametrization of the PDFs together with the values nic = 



1.3 GeV and A = 326 MeV [Q.^ Unless otherwise stated, we use n = sj 4m^ + throughout 
this paper. 

For comparison, the panel [d) of Fig. |3| shows the ^ dependence of the QCD coiTcction 
factor for the transverse structure function, K(x^Q^) = Fj^^ /Fj*^. One can see that sizable radia- 
tive corrections to the structure functions Ft{x,Q^) and Fl{x,Q^) cancel each other in their ratio 

Q^) = Fl/Ft with good accuracy. As a result, the NLO contributions to the ratio R{x, Q^) are 
less than 10% for x > 10""*. 

Another remarkable property of the Callan-Gross ratio closely related to fast perturbative con- 
vergence is its parametric stability.^ Our analysis shows that the fixed-order predictions for the 
ratio R{x, Q^) are less sensitive to standard uncertainties in the QCD input parameters than the cor- 
responding ones for the production cross sections. For instance, sufficiently above the production 
threshold, changes of jU in the range {l/2)^/4m^ + < ^ < 1\J Am}. + only lead to 10% vari- 
ations of R{x,Q^) at NLO. For comparison, at x = 0.1 and = 4.4, such changes of jX affect the 
NLO predictions for the quantities Ft{x, Q^) and R{x, Q^) in charm leptoproduction by more than 
100% and less than 10%, respectively. 



^Note that we convolve the NLO CTEQ6M distribution functions with both the LO and NLO partonic cross sections 
that makes it possible to estimate directly the degree of stability of the pQCD predictions under radiative corrections. 

^Of course, parametric stability of the fixed-order results does not imply a fast convergence of the corresponding 
series. However, a fast convergent series must be parametrically stable. In particular, it must exhibit feeble jXf and jXR 
dependences. 
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Figure 3: (a), {b) and (c) panels: dependence of the LO (solid curves) and NLO (dashed curves) 
predictions for the Callan-Gross ratio, R{x,Q^) = Fi/Ft, in charm leptoproduction at x = 10^', 10^^ and 
10^^. {d) panel: dependence of the /iT factor for the transverse structure function, K[x, Q^) = Fj^'^/Fj'^, 
at the same values of x. 



Keeping the value of the variable fixed, we analyze the dependence of the pQCD pre- 
dictions on the uncertainties in the heavy-quark mass. We observe that changes of the charm- 
quark mass in the interval 1.3 GeV < nic < 1.7 GeV affect the Callan-Gross ratio by 2%-3% at 
= 10 GeV^ and x < 10^ ^ The corresponding variations of the structure functions Ft{x,Q^) 



and Fl{x,Q^) are about 20%. We also verify that the recent CTEQ versions [38, 39] of the PDFs 
lead to NLO predictions for R{x,Q^) that coincide with each other with an accuracy of about 5% 

at 10"^ <x< m-\ 



3. Soft-gluon corrections to the azimuthal asymmetry A (x, Q^) at NLO 

Presently, the exact NLO predictions for the azimuth dependent structure function Fa{x,Q^) 
are not available. For this reason, we consider the NLO predictions for the azimuthal cos(2(p) 
asymmetry within the soft-gluon approximation. For the reader's convenience, we collect the final 
results for the parton-level GF cross sections to the next-to-leading logarithmic (NLL) accuracy. 
More details may be found in Refs. [0 [l|, 0, ||]. 

At NLO, photon-gluon fusion receives contributions from the virtual ^(aemO;^^) coiTcctions to 



the Bom process (2.1 1 and from real-gluon emission, 



f {q) + sikg) -^Q{pq) + Q{pq) + g{pg 



(3.1) 
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The partonic invariants describing the single-particle inclusive (IPI) kinematics are 



s 



2q-kg = s + Q 

s' + ti+Ui, 



u\ = [q-PQ? 



- m 

■ 17?' 



(3.2) 



where C, is defined through kg = C,p and ^'4 measures the inelasticity of the reaction (3.1). The 
corresponding IPI hadron-level variables describing the reaction ( |1.1[ ) are 



S' = 2q-p = S+Q\ 
S4 = S' + Ti+Ui, 



Ti 



(p-Pq) -m^ 

{q-Paf-m^ 



(3.3) 



The exact NLO calculations of unpolarized heavy-quark production [[T]-^] show that, near 
the partonic threshold, a strong logarithmic enhancement of the cross sections takes place in the 
collineai; \pg.T \ — )■ 0, and soft, \pg\ — > 0, limits. This threshold (or soft-gluon) enhancement is of 
universal nature in perturbation theory and originates from an incomplete cancellation of the soft 
and collinear singularities between the loop and the bremsstrahlung contributions. Large leading 
and next-to-leading threshold logarithms can be resummed to all orders of the perturbative expan- 



sion using the appropriate evolution equations [|40||. The analytic results for the resummed cross 
sections are ill-defined due to the Landau pole in the coupling constant a^. However, if one con- 
siders the obtained expressions as generating functional and re-expands them at fixed order in a^, 
no divergences associated with the Landau pole are encountered. 



Soft-gluon resummation for the photon-gluon fusion was performed in Ref. [ ]1Q| ] and confirmed 
in Refs. [jTs], |l^. To NLL accuracy, the perturbative expansion for the partonic cross sections, 
d?ak{s' ,t\,u\) / [dty Au\) (k = T,L,A,I), can be written in factorized form as 



dtidui 



{s\tuu,)=Bl°^{s'A,u, 



5{s +ti+Ui) + Yj 



n=\ 



71 



K^''\s\tuu,) 



(3.4) 



The functions K^''\s' ,?i,Mi) in Eq. ( p^ originate from the collinear and soft limits. Since 
the azimuthal angle <p is the same for both Y'g and QQ center-of-mass systems in these limits, 
the functions K^'^\s\t\^ui) are also the same for all du, k = T,L,A,I. At NLO, the soft-gluon 



coixections to NLL accuracy in the MS scheme read [10] 



In [s^/m^ 



J + 



1+ln 



h 



+ 5(^4 In— r In 



(3.5) 



In Eq. (|3|), Ca = Nc, Cp = {N^ - \)/{2Nc), Nc is the number of quark colors, and = (1 - 
2m^/s){\n[{\ -j8J/(l +j3-)]+i7r} with j3- = ^Jl-Am^/s. The single-particle inclusive "plus" 
distributions are defined by 



In' {s4/m^) 
sa 



lim 



In' (54 / 



-e{sA-e) + 



' -In'+'A 



/ + 1 



5(^4) 



(3.6) 
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For any sufficiently regular test function h{si\), Eq. implies that 









-L 

+ 



' d.4 [h{s,) - h{0)] ^ihlnfl + J_/,(o) in'+i 



1+1 



Jo 

(3.7) 

Standard NLL soft-gluon approximation allows us to determine unambiguously only the sin- 



gular behavior of the cross sections defined by Eq. (3.6). To fix the dependence of the Bom- 
level distributions Bf°^{s' ,ti,U[)\^^ -s -s'-t ^1" QA^^ ^^'^ method proposed in and 



based on comparison of the soft-gluon predictions with the exact NLO results. According to [|20||, 



X4 



Ui 



s' + h 



1 



^4 



(3.8) 



where the LO GF differential distributions ^^&df'^{s\t\)/dt\ are 



Ke^Qa^^a, i — + — +4 



2 / 

s m s 



t\U\ 



tiU\ s' 



2 82" 



5 m^s 



t\U\ 



2 / 

2 . I s m s 

KegasmOCs 4 -7 - - — 



2 ^ fT^ fhu-is 

7r4aema.4V2^ ( — ^ 



■ m 



«1 =— i'— /l 

— + ^ 

2 \ ^1 ~ "1 



(3.9) 



22^ 2m2V 



1 



tlU\ 



U{= — S' — t{ 



Comparison with the exact NLO results given by Eqs. (4.7) and (4.8) in Ref. [y] indicates 
that the usage of the distributions Bf°'^{s' ,t\,ui) defined by Eqs. (3^) and (3^) in present pa- 



per provides an accurate account of the logarithmic contributions originating from collinear gluon 
emission. Numerical analysis shows that Eqs. (3.8) and ( |3.9D render it possible to describe with 
good accuracy the exact NLO predictions for the functions dj\s') and di\s') near the threshold 
at relatively low virtualities ~ [20|.^ 

Our results for the x distribution of the azimuthal cos(2<p) asymmetry, A{x, Q^) = IxFa/Fi, in 
chaim leptoproduction at fixed values of are presented in the left panel of Fig. Q For comparison, 
the K factor, K{x, Q^) = F2^'^ /F]^'^, for the structure function F2 at the same values of t, is shown 
in the right panel of Fig. ^ One can see that the sizable soft-gluon corrections to the production 
cross sections affect the Bom predictions for A(.x:, Q^) at NLO very little, by a few percent only. 

4. Analytic LO results for R{x, Q^) and A(;c, Q^) at low x 

Since the ratios R{x,Q^) and A{x,Q^) are perturbatively stable, it makes sense to provide 
the LO hadron-level predictions for these quantities in analytic form. In this Section, we derive 



^Note that soft-gluon approximation is unreliable for high 2> m^. 
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Figure 4: Left panel: LO (solid lines) and NLO (dashed lines) soft-gluon predictions for the x dependence 
of the azimuthal cos(2^) asymmetry, A{x,Q^) = 2xFa/F2, in charm leptoproduction at = 1 and 5. Right 
panel: x dependence of the K factor, K[x, Q^) = Fj'^'^^/Fj^^, at the same values of E, . 



compact \ow-x approximation formulae for the azimuthal cos(2<p) asymmetry and the quantity 
/?2(-^5 Q^) closely related to the Callan-Gross ratio R{x, Q^), 



Ri(x,Q')=1x%x,Q^) 



R{x,Q^ 



(4.1) 



l+R{x,Q^)' 

We will see below that our obtained results may be useful in the extraction of the structure functions 
(k = 2,L,A,I) from experimentally measurable reduced cross sections. 
To obtain the hadron-level predictions, we convolute the LO partonic cross sections given by 
Eqs. (2.2) with the low-.x asymptotics of the gluon PDF: 

1 



A+S ' 



(4.2) 



The value of 5 in Eq. ( p^^ ) is a matter of discussion. The simplest choice, 5 = 0, leads to 
a non-singular behavior of the structure functions for x ^ 0. Another extreme value, 5 = 1/2, 
historically originates from the BFKL resummation of the leading powers of ln(l/.x) [^-|^]. In 
reality, 5 is a function of (for an experimental review, see Refs. [|T| Hi). Theoretically, the 
dependence of 5 is calculated using the DGLAP evolution equations 



32]. 



First, we derive an analytic low-.x formula for the ratio /?2 (2^) = ^2 (-^ ~^ 0)2^) with ar- 
bitrary values of 5 in terms of the Gauss hypergeometric function. Our result has the following 
form: 



1+4A. 



1 I g(i-g^) 

^ (2+5)(3+5) 



1+4A. 



(l+4A)(4-5- 



10 

3+5 



(4.3) 



1+4A. 



where A is defined in Eq. ( p3| ) and the function (£>{r,z) is 



<I>(r,z) 



r(i/2)r(i + r) 



1 



2/^1(2^1 + ^;- 



■r;z 



l+r r(3/2 + r) 

The hypergeometric function 2Fi(a,b;c;z) has the following series expansion: 
r7 f h ^ r(c) f ria + n)r{b + n) z'^ 



(4.4) 



(4.5) 
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Figure 5: LO low-x predictions for the ratio R2{x,Q^) = 2xFi/F2 in charm leptoproduction. Left panel: 
Asymptotic ratios ^2"^ {Q^) (gray points) and /Jj'^^' (2^) (black points), as well as CTEQ6L predictions for 
Riix.Q^) at jc = 10"^ 10"^ and lO"'^. Right panel: Asymptotic ratio r[^\q^) at 5 = 0, 0.2, 0.3, 0.4 and 
0.5. 



In Fig. ^, we investigate the obtained result (4.3) for 7?2^'(2^). The left panel of Fig. || shows 
the ratio /?2^- (2^) as functions of for two extreme cases, 5 = and 1/2. One can see that the 
difference between these quantities varies slowly from 20% at low to 10% at high Q^. For 
comparison, also the LO results for R2{x,Q^) calculated at several values of x using the CTEQ6L 
gluon PDF p8[ ] ai^e shown. We observe that, for x — )• 0, the CTEQ6L predictions converge to the 
function R2 {Q^) practically in the entire region of Q^. We have verified that the similar situation 
takes also place for other recent CTEQ PDF versions [38, 39 1. In the right panel of Fig. B| the d 
dependence of the asymptotic ratio /?2^'(2^) is investigated. One can see that the ratio /?2^'(2^) 
rapidly converges to the function R2^^\q'^) for 5 > 0.2. In particular, the relative difference 
between R2'^\q^) and R2'^\q^) varies slowly from 6% at low to 2% at high Q^. 

As mentioned above, the dependence of the parameter 8 is determined with the help of the 
DGLAP evolution. However, our analysis shows that hadron-level predictions for /?2^^(x — )• 0, Q^) 
depend weakly on 5 practically in the entire region of for 0.2 < 5 < 0.9. For this reason, it 
makes sense to consider the ratio ^2^' (2^) in particular case of 5 = 1 /2. The result is: 

/?(l/2).^2^ ^ [3+4A(13 + 32A)]£(l/(l+4A))-4A {9 + 32X)K{\/{\+AX)) 

2 ' 1+4A (-37 + 72A)£(l/(l+4A)) + 2(23-36A)^(l/(l+4A)) ' 

where the functions K{y) and E{y) are the complete elliptic integrals of the first and second kinds 
defined as 

K{y)=f^=== f E{y)= fdtJ^—^. (4.7) 

^0 ^[\-fl)[l-yfi) Jo V 



One can see from Fig. g that our simple formula ( p^ with 5 = 1/2 (i.e., without any evolution) 
describes with good accuracy the low-x CTEQ results for R2{x, Q^). We conclude that the hadron- 
level predictions for R2{x ^ 0,Q^) are stable not only under the NLO corrections to the partonic 
cross sections, but also under the DGLAP evolution of the gluon PDF. 

Then we calculate and investigate the LO hadron-level predictions for the azimuthal cos(2(p) 
asymmetry in the limit of x — ^ 0. Our result for the quantity A^^^ {Q^) = A^^) (x — 0, 2^) has the 
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1 10 100 1000 10000 1 10 100 1000 10000 

if = Q2 / if = Q2 / 



Figure 6: LO low-x predictions for the ratio A(x,(2^) = 2xFa/F2 in charm leptoproduction. Left panel: 
Asymptotic ratios A'"' {Q^) (gray points) and A^'/^^ {Q^) (black points), as well as CTEQ6L predictions for 
A{x,Q^) atx= 10"^ 10"^ and 10"'*. Right panel: Asymptotic ratio A (^^(e^) at 5 =0, 0.2, 0.3,0.4 and 0.5. 



following form: 



1 + 



(2+6-)(3+5) 



<J> (5, T^) - (1 +4A) (4 - 5 - 3!^^) O (1 + 5, ^) 



(4.8) 



Our analysis presented in Fig. ^ shows that the quantity A'^^\Q^) defined by Eq. ( |4.8| ) has the 
properties very similar to the ones demonstrated by the ratio R^\q^). In particular, one can see 
from Fig. ^ that the hadron-level predictions for A^^\Q^) depend weakly on 5 practically in the 
entire region of for 5 > 0.2. So, the azimuthal cos(2<p) asymmetry A{x — )• 0, Q^) is also stable 
under the DGLAP evolution of the gluon PDF. 

Let us now discuss how the obtained analytic results may be used in the extraction of the 
structure functions (k = 2,L,A,I) from experimentally measurable quantities. Usually, it is the 
so-called "reduced cross section", a(x, Q^), that can directly be measured in DIS experiments: 



1 



l + {l-yy2nai^dxdQ^ 



2xy^ 



i + {i-yr 

„2 



1 



y 



(4.9) 
(4.10) 



In earlier HERA analyses of charm and bottom electroproduction, the corresponding longitudinal 
structure functions were taken to be zero for simplicity. In this case, a(x, q2) =F2{x,Q^). In recent 
papers [|^, 0], the structure function F2(x, 2^) is evaluated from the reduced cross section (1.9) 
where the longitudinal structure function Fi{x,Q^) is estimated from the NLO QCD expectations. 



Instead of this rather cumbersome procedure, we propose to use the expression ( [1.10[ ) with the 
quantity R2{x, Q^) defined by the analytic LO expressions (4.3) or ( 1.6). This simplifies the extrac- 



tion of F2{x,Q^) from measurements of d{x,Q^) but does not affect the accuracy of the result in 
practice because of perturbative stability of the ratio R2{x, Q^). 

In Ref. [pQ[], we used the analytic expressions ( |4.3[ ) and (4.6) for the extraction of the structure 
functions F2{x,Q^) and F2{x,Q^) from the HERA measurements of the reduced cross sections 



a''(x, Q ) and a (x, Q ), respectively. It was demonstrated that our LO fomiula (1.6) for R2{x, Q 
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with 8 
laboration [I 



1/2 usefully reproduces the results for F2{x, Q^) and Q^) obtained by the HI Col- 
3^] with the help of the NLO evaluation of Fl{x,Q^) 



In particular, the results of 
our analysis of the HERA data on the chami electroproduction are collected in Table |l|. In our cal- 
culations, the value = 1.3 GeV for the charm quark mass is used. The LO predictions, F2(L0), 
for the case of 5 = 0.5 are presented and compared with the NLO values, f2(NL0), obtained in 
the HI analysis |3^. One can see that our LO predictions agree with the NLO results with an 



accuracy better than 1%. 



Table 1: Values of (jCjQ^) extracted from the HERA measurements of a'^{x,Q^) at low [34| and high 
[ p3| | (in GeV^) for various values of x (in units of 10^^). The NLO HI results p3| , p4[ | are compared with 
the LO predictions corresponding to the case of 5 = 0.5. 





X 


y 


a' 


Error 


F2'^(NL0) 




(GeV2) 


(xlO-3) 






(%) 


HI 


5 = 0.5 


12 


0.197 


0.600 


0.412 


18 


0.435 ±0.078 


0.435 ±0.078 


12 


0.800 


0.148 


0.185 


13 


0.186 ±0.024 


0.1 85 ±0.024 


25 


0.500 


0.492 


0.318 


13 


0.331 ±0.043 


0.331 ±0.043 


25 


2.000 


0.123 


0.212 


10 


0.212 ±0.021 


0.212 ±0.021 


60 


2.000 


0.295 


0.364 


10 


0.369 ±0.040 


0.369 ±0.040 


60 


5.000 


0.118 


0.200 


12 


0.201 ±0.024 


0.200 ±0.024 


200 


0.500 


0.394 


0.197 


23 


0.202 ±0.046 


0.202 ±0.046 


200 


1.300 


0.151 


0.130 


24 


0.131 ±0.032 


0.130±0.031 


650 


1.300 


0.492 


0.206 


27 


0.213 ±0.057 


0.213 ±0.057 


650 


3.200 


0.200 


0.091 


31 


0.092 ±0.028 


0.091 ±0.028 



The structure functions Fa and Fj can be extracted from the <p-dependent DIS cross section, 
d^GiN 2al„ r 1 



dxdQ^dq) l-e 



2x 



F2{x,Q')-{1-£)Fl{x,Q') 



(4.11) 



where e 



+eFAix,Q'^)cos2(p + 2y/£{l + e)Fi{x,Q'^)cos(p , 
For this purpose, one should measure the first moments of the cos((p) and 



cos(2<p) distributions defined as 

J^-dcpcosncp^^ix,Q\cp) 
{cosn(p}{x,Q^^ — ^ 



(4.12) 



Using Eq. (4.11), we obtain: 



{cos2(p){x,Q^) 
{cos(p){x,Q^) 



1 £A{x,Q^ 



Fa 



2 1-(l-£)/?2(x,e2) 
^e{\ + e)A,{x,Q^) 
l-{l-e)R2{x,Q^)' 



A{x,Q^) = 2x-^{x,Q^), 
A,ix,Q^) = 2x^{x,Q^). 



(4.13) 
(4.14) 
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One can see from Eqs. (4.13) and (4.14) that, using the perturbatively stable predictions (4.3) 



for R2{x,Q ), we will be able to determine the structure functions Fa{x,Q ) and Fi{x,Q ) from 



future data on the moments (cos2<p) and (cos<p). On the other hand, according to Eq. ( |4.13| ), 
the analytic results (^) and ( |4.8[ ) for the quantities R2{x,Q^) and A{x,Q^) provide us with the 
perturbatively stable predictions for (cos2<p) which may be directly tested in experiment. 

So, our obtained analytic and perturbatively stable predictions for the ratios R2ix,Q^) and 
A{x,Q^) will simplify both the extraction of structure functions from the measurable (p-dependent 
cross section ( [4.1 1[ ) and the test of self-consistency of the extraction procedure. 

5. Conclusion 

We conclude by summarizing our main observations. In the present paper, we studied the ra- 
diative coixections to the Callan-Gross ratio, R{x, Q^), and azimuthal cos(2(p) asymmetry, A{x^ Q^), 
in heavy-quark leptoproduction. It turned out that large (especially, at non-small x) radiative cor- 
rections to the structure functions cancel each other in their ratios R{x, Q^) = Fi/Ft and Q^) = 
2xFa/Fa with good accuracy. As a result, the NLO contributions to the ratios R{x, Q^) and Q^) 
are less than 10% in a wide region of the variables x and Q^. Our analysis shows that, sufficiently 
above the production threshold, the pQCD predictions for R{x,Q^) and A{x,Q^) are insensitive 
(to within ten percent) to standard uncertainties in the QCD input parameters and to the DGLAP 
evolution of PDFs. We conclude that, unlike the production cross sections, the Callan-Gross ratio 
and cos(2(p) asymmetry in heavy-quark leptoproduction are quantitatively well defined in pQCD. 
Measurements of the quantities R{x, Q^) wAA{x, Q^) in charm and bottom leptoproduction would 
provide a good test of the conventional parton model based on pQCD. 

As to the experimental aspects, we propose to exploit the observed perturbative stability of 
the Callan-Gross ratio and azimuthal asymmetry in the extraction of the structure functions from 
the experimentally measurable reduced cross sections. For this purpose, we provided compact LO 
hadron-level formulae for the ratios R2{x, Q^) = 2xFi/F2 = R/{l+R) and A{x,Q^) = 2xFa/F2 in 
the limit x^O. We demonstrated that these analytic expressions usefully reproduce the results for 
F2{x, Q^) and F2{x, Q^) obtained by the HI Collaboration [33, |^ with the help of the more cum- 
bersome NLO evaluation of Fi{x,Q^). Our obtained predictions will also be useful in extraction 
of the azimuthal asymmetries from the incoming COMPASS results as well as from future data on 



heavy-quark leptoproduction at the proposed EIC [ [35[ j and LHeC [ [36[ ] colliders at BNL/JLab and 
CERN, coiTcspondingly. 
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